Bulk and shear viscosities of matter created in relativistic heavy-ion collisions 
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We study the effects of shear and bulk viscosities in the hadronic phase on the expansion of the 
fireball and on the particle production in relativistic heavy ion collisions. Comparing simulation 
with or without viscosity in the hadronic matter we find that elliptic flow observables strongly 
dependent on dissipative effects in the late stage. On the other hand, interferometry radii are 
sensitive, through the early transverse flow, on the value of the viscosity at high temperatures. We 
present first calculations including the effects of bulk viscosity in the hadronic phase and in the 
hadron emission. We find them important in obtaining a small freeze-out temperature consistent 
with the measured transverse momentum spectra and elliptic flow of identified particles. 
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I. INTRODUCTION 

The matter created in relativistic heavy-ion collisions 
at the BNL Relativistic Heavy Ion Collider (RHIC) is a 
dense strongly interacting fluid [1-4] . The observation of 
strong collective transverse and elliptic flows is an indi- 
cation that the system behaves as a fluid. To model the 
dynamics of such a system relativistic hydrodynamics of 
a perfect fluid has been successfully applied [5-11]. The 
fireball expands and cools down, until some freeze-out 
temperature is reached, after which particles are emitted 
from a freeze-out hypersurface. Final particle spectra, 
to be compared with experimental data, are obtained af- 
ter resonance decays. Transverse momentum spectra in 
the azimuthal angle at central rapidity are written as an 
expansion in Fourier coefficients 
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(1.1) 



The form of the observed transverse momentum spec- 
and the elliptic flow coefficient vi can be 



dN 



2irp±dp±dy 

described using a convolution of the thermal emission 
of particles with the collective velocity of the fluid itself 
[12, 13]. 

Due to the rapid expansion of the dense system cre- 
ated in relativistic heavy-ion collisions deviations from 
local equilibrium can be important. For the modeling of 
the expansion of the fireball it means that viscous rela- 
tivistic hydrodynamics should be used [14-19]. A con- 
sistent causal scheme requires the use of second order 
viscous equations [20] . Most of the existing applications 
of viscous hydrodynamics in heavy-ion collisions consider 
shear viscosity only. The value of the ratio r)/s of the 
shear viscosity coefficient to the entropy density is an im- 
portant characteristic of the strongly interacting medium 
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created in the course of the collision [21, 22]. The extrac- 
tion of the shear viscosity coefficient from the measured 
elliptic flow could give valuable information [17, 23]. The 
main source of uncertainty in the analysis lies in the as- 
sumption of the initial eccentricity of the source at a 
given impact parameter [23]. The causes influencing the 
initial shape of the source include : different underlying 
models of the initial density, Color Glass Condensate or 
Glauber Model [24] , the inclusion of binary collision con- 
tributions [25], possible initial fluctuations of the shape 
(standard versus participant eccentricity) [26, 27] or a 
core-corona effect, where only the dense part of the source 
evolves collectively [28] . The elliptic flow of the bulk of 
the matter is generated in the early stages of the colli- 
sion. However, the final elliptic flow of observed hadrons 
is modified in the hadron gas phase of the expansion, due 
to rescattering and resonance decays [10, 29-31]. This 
is true both for the elliptic flow of charged particles as 
well as of identified particles [10, 32]. In particular, to 
reproduce the observed splitting between pions and pro- 
tons in the transverse momentum dependence of V2 a late 
freeze-out or a hadronic cascade stage are required in the 
evolution. 

The role of dissipation in the hadronic phase must be 
assessed before a reliable estimate of viscosity in the 
(quark-gluon plasma) QGP phase can be made. Al- 
though the importance of the difference of the viscosity 
coefficients in the hadronic and plasma phases has been 
discussed [29], most of the existing hydrodynamic simu- 
lations applied to heavy ion collisions use a constant r\j s 
ratio through the evolution. In this paper we study the 
effect of viscosity in the hydrodynamic evolution below 
the transition temperature on the final elliptic flow, spec- 
tra and Hanburry Brown- Twiss (HBT) correlation radii. 
In particular we analyze the observable differences in soft 
momenta observables between two extreme assumptions 
on the shear viscosity in the plasma phase {i]/s = or 
0.16), after hadronic dissipation is taken into account. 
We show that the effect of dissipation in the hadronic 
phase strongly reduces the sensitivity of the elliptic flow 
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measure on the value of the viscosity in the early QGP 
phase of the expansion. If the hadronic phase in the 
expansion is dilute enough, a cascade after burner can 
used after an early freeze-out of the fluid [31, 33-35]. 
Alternatively a longer hydrodynamic evolution can be 
used with a hadronic equation of state below the transi- 
tion temperature. This paper studies the effect of such 
a longer hydrodynamic evolution in the hadronic phase 
using viscous hydrodynamics. We use a moderate value 
of rj/s = 0.1 in the hadronic phase and a bulk viscosity 
Q/s = 0.03-0.04. We show that even such small values 
of viscosities in the hydrodynamic evolution in the late 
phase of the collision are important for the final elliptic 
flow, and that with such assumptions we can reproduce 
many experimental observations. We calculate also the 
HBT radii after a hydrodynamic evolution with different 
viscosities in the QGP and hadron gas phases. 



II. SHEAR AND BULK VISCOSITIES 

Besides the ideal fluid expansion we consider three 
other scenarios for the shear and bulk viscosities in the 
hot matter. The general idea is that the shear viscosity 
in the hadronic and QGP phases could be very different. 
Moreover if the shear viscosity in the hadronic phase is 
non-zero, it could be accompanied by non-negligible bulk 
viscosity. The formula for the temperature dependence 
of the ratio of the shear viscosity to the entropy is taken 
in the form 



,(T)^£f HG (T) + (l-f HG (T))^- (2.1) 



where the function /hg(^) = 

l/(exp((T-T c )/AT) + l) cuts-off the hadron gas 
viscosity above T c = 170MeV ( AT = 8MeV ). 
fio W {T) = 1/ (exp ((T low - T)/AT) + 1) is introduced to 
cut-off viscosity effects below Ti ow = 80MeV to improve 
numerical stability. Depending on the chosen values of 
the viscosities in the hadronic matter and in QGP we 
consider four different scenarios (Table I). The tem- 
perature dependence of the viscosities corresponding to 
viscous scenarios in the Table are shown in Fig. 1. The 
scenario denoted as vHG assumes that only viscosity 
in the hadronic phase is non-zero. vQGP is taken for 
illustration to show how the dissipative phenomena in 
the plasma alone could influence the final observables. 
The scenario vQGP+vHG is the most general (with a 
suitable choice of t}hg and t]qgp)- The two scenarios 
vHG and vQGH+vHG differ by the choice of the shear 
viscosity coefficient in the plasma. The comparison of 
these two different scenarios is one of the motivations of 
this investigation, namely to test how sensitive the final 
observables are to the assumed viscosity in the plasma 
phase (from rj/s = to rj/s = 0.16), when dissipation 
in the hadronic phase occurs afterwards. These two 
scenarios scenarios represent two extreme assumptions 
on the temperature dependence of the ratio rj/s, i.e. 
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VQGP 


T F (MeV) 
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140 


vHG 


0.1 





150 


vQGP 





0.16 


130 


vHG+vQGP 


0.1 


0.16 


135 



TABLE I: Viscosity parameters used in the four calculations 
presented in the paper. The last column contains the freeze- 
out temperature that reproduces best pion spectra in each 
case. 



increasing or decreasing when switching from the QGP 
to hadronic matter. Existing viscous hydrodynamic 
simulations assume a constant rj/s as function of temper- 
ature, microscopic estimates suggest a (strong) increase 
of rj/s when decreasing the temperature. We test a 
scenario with a moderate increase of rj/s at T c and also 
another extreme scenario where the reverse happens and 
show that the results are in fact very similar and that 
reproducing experimental data requires a small value of 
rjHG for any QGP viscosity. 

The near equilibrium regime in a dilute gas of interact- 
ing hadrons can be described using the Boltzmann equa- 
tion. Estimates of shear viscosity with hadronic cross 
sections or chiral models lead to a large value rj/s ~ 1 
for temperatures T ~ 150MeV [36-40]. The large value 
rj/s ~ 1 in the hadronic phase seems to contradict exist- 
ing fits of the data using viscous hydrodynamics, where 
rj/s = 0.08-0.16 depending on the initial eccentricity [17]. 
Also such a large value of the viscosity coefficient would 
simply mean that the viscous hydrodynamics cannot be 
applied. Shear viscosity could be significantly reduced 
if the number of hadronic states increases near T c [41]. 
At temperatures close to the transition temperature the 
description of the dense medium using a transport equa- 
tion approach involving hadrons with vacuum properties 
could break down. On the other hand, in microscopic 
models the bulk viscosity is estimated to be much smaller 
C/s ~ 0.03-0.05 [41, 42]. 

In this paper, we use relativistic viscous hydrodynam- 
ics to model the dense hot medium on the low tempera- 
ture side of the transition temperature. The equation of 
state of matter for T < T c is approximated as the hadron 
gas equation of state involving 371 known hadrons and 
resonances [43] . This equation of state can be smoothly 
connected to the equation of state calculated in lattice 
QCD at higher temperatures. The final equation of state 
leads to a correct description of spectra and HBT radii 
in ideal fluid hydrodynamics [11, 44]. Shear viscosity 
in the hadronic phase is treated as a free parameter in 
our calculation. A simple estimate of the viscosities can 
be obtained in the relaxation time approximation [45- 
47] . Starting from the Boltzmann equation for the phase 
space distribution distribution /„ of particle species n 
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where Sf n = f n — is the deviation from the equilib- 
rium (Bose-Einstein or Fermi-Dirac) distribution f® — 
cx P (p>'n^( a: )/T)±i ' u n( x ) is tne local nuid four- velocity, 
and thg? is the relaxation time (the same for all particle 
species). In the local rest frame we have 

Sfn = Y§f n (1 ± O (PVO^ - c 2 s E 2 d^) . (2.3) 

Calculating the stress corrections to the energy momen- 
tum tensor 

ST MV = ^ + nA ^ = J2 J J^z^Sfn > (2-4) 

where A Ml/ = g^ v — u^u u , we get for the stress tensor in 
the local rest frame 

« thg \- /" d 3 p p l p>p k p l „ „ „ 

n ~—l^J ( 2 ^)3 g2 A. I 1 ± /n) g ( 2 -5) 

and 

(2.6) 

with 

CT Q £ = i ^V a ll(j + V^U Q - ^AapdfiU'^J . (2.7) 

Comparing with the first order expressions for the stress 
tensor n^u = 2^^, II = —^d^u^ we have 

and 

.(2.9) 

In the modelling of heavy ion collisions we are inter- 
ested in the properties of the hadronic matter in a tem- 
perature range from the freeze-out temperature Tp > 
130MeV to the transition temperature T c = 170MeV. 
Performing the sums over the resonances used in the 
calculation of the hadronic matter equation of state, 
i.e. the resonances listed in the SHARE program [48], 
we can relate the viscosity coefficient to the relaxation 
time. Assuming a constant shear viscosity to entropy ra- 
tio rj/s = 0.1 between 80 and 170MeV, the relaxation 
time thg changes weakly in the range 0.8-1.2fm/c. We 
could assume instead a different dependence of rj/s or of 
thg on the temperature, but these details do not mat- 
ter much. It turns out that it is the value at freeze-out 
that is the most important. Our choice corresponds to 
thg — lfm/c at T = 150MeV, but other values of the pa- 
rameters could be tested in more extensive sets of model 
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FIG. 1: (Color online) Temperature dependence of the ratio 
of shear and bulk viscosities to the entropy. The solid, dotted 
and dash-dotted lines represent the shear viscosity for the 
vHG, vQGP and vQGP+vHG scenarios. The dashed line 
represent the bulk viscosity, effective only in the vHG and 
vQGP+vHG scenarios. 

calculations. The assumed small hadronic shear viscos- 
ity is not motivated by microscopic model estimates, that 
would suggest a larger value. It is rather motivated by 
existing viscous hydrodynamic calculations [17, 23], in- 
dicating that the average viscosity in the hadronic and 
QGP phases is small. 

The bulk viscosity of the hadronic matter at high den- 
sity is another not very well controlled parameter. It is 
expected that in the deconhned phase the bulk viscosity 
coefficient is negligible. On the other hand, a sharp rise 
of the bulk viscosity has been predicted [49] around T c . 
If the effect of the bulk viscosity at T c is large the flow 
could be modified [50, 51] or could even become unstable 
leading to the fragmentation of the fireball [52] . On the 
other hand, the rise of the bulk viscosity near T c could 
be accompanied by critical slowing down, which leads 
to an increase of the dynamical bulk viscosity relaxation 
time Tn, delaying the onset and effectively diminishing 
bulk viscosity effects. 2 + ID hydrodynamic simulations 
indicate that by the time the expanding system reaches 
T c substantial amount of transverse flow has already set 
in [50] and the effects of the rising bulk viscosity at the 
critical temperature is reduced and the agreement of the 
calculation with the data is not spoiled. In this paper we 
do not take into account bulk viscosity near the phase 
transition. 

In the hadron gas phase the bulk viscosity could be 
quite substantial, as particle masses get comparable to 
the temperature. Bulk viscosity can be estimated in the 
relaxation time approximation from Eq. (2.9). The re- 
sulting C/s corresponding to t]hg/ s — 0.1 is shown in 
Fig. 1. Since the relaxation time formulas use physical 
hadrons, we restrict the temperature range for the cal- 
culation in Eqs. 2.8 and 2.9 to the hadronic phase, tak- 
ing 2(T) = f low (T)^f HG (T). The shear viscosity at 
larger temperatures (1 — fHG(T)) VQ ° P is not generated 
through hadronic processes. In this paper we are inter- 
ested in the effects of bulk viscosity in the late stages, 
and therefore we do not take into account possible bulk 
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viscosity of non-hadronic origin. At temperatures around 
150MeV we have C/ s — 0.035. For our estimate of vis- 
cosities using a relaxation time formula with thg of the 
order of lfm/c we obtain the bulk viscosity similar as in 
microscopic models [42] but the shear viscosity is signif- 
icantly smaller than in most estimates [36-40] . To check 
this assumption we performed also a calculation with the 
same bulk viscosity but increasing t]hg / s to 0.24. This 
would mean that relaxation time formulas do not apply. 
We find that the assumed value of the shear viscosity 
Vhg/ s — 0.24, which is still smaller than microscopic 
estimates, gives already a too strong suppression of the 
elliptic flow. 

Eq. (2.2) defines nonequilibrium corrections to the dis- 
tribution function. The corrections from bulk viscosity 
cannot be taken in the form of the Grad's expansion [53] . 
From Eq. (2.2) we get for the corrections from bulk vis- 
cosity IT [47] 



5f b n ulk = C bulk f° n (1 ± /°) ( c\E 

in the local rest frame, with 
1 1 x - f d 3 p m 2 
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(2.10) 
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(2.11) 

The deviation from equilibrium due to the stress correc- 
tions from shear viscosity are taken in the form [15, 54] 



S j 'shear - /„ (l ± /„) 



1 



2T 2 (e+p) 



(2.12) 



with e the local energy density and p the pressure. It 
must be noted that more general forms of the nonequi- 
librium corrections are possible for multicomponent sys- 
tems or for species dependent relaxation times [55] . 



III. VISCOUS HYDRODYNAMIC EVOLUTION 
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FIG. 2: (Color online) Time dependence of the entropy scaled 
by the initial entropy in the ideal fluid calculation. The solid, 
dashed, dash-dotted and dashed lines represent results from 
the ideal fluid, vQGP, vQGP+vHG and vHG calculations re- 
spectively. 
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FIG. 3: (Color online) Freeze-out hypersurfaces T(t,x,y = 
0) = Tf at impact parameter b — 2.2fm. The solid, dashed, 
dash-dotted and dashed lines represent hypersurfaces in the 
ideal fluid, vQGP, vQGP+vHG and vHG calculations respec- 
tively. 



The hydrodynamic equations 



and 



dyT^ = 



(3.1) 



are solved in 2 + ldimensions, assuming boost invariance 
in the longitudinal direction. The energy momentum ten- 
sor 

T» v = (e + p)u^u" - pg» v + tt^ + IW (3.2) 

is composed of the ideal fluid part and stress shear and 
bulk viscosity corrections tt and II. The viscous correc- 
tions in the second order Israel-Steward viscous hydro- 
dynamics are solutions of the dynamical equations [20] 



(3.3) 



u 7 a 7 n = 



-C<9 7 u 7 -n 1 rr 



i rn 



(3.4) 



and assume 
= 7r»»(7a) = 



We take for the relaxation time t w = 
T~n =T W . The initial conditions are tt xx (to) 
gL, it x v(to) = and n(r ) = 0. The details of the 
choice of initial conditions and ttj are not crucial, as the 
bulk viscosity correction gets rapidly close to the Navier- 
Stokes value, and anyway its influence on the dynamics 
itself is small. 

For the energy density profile in the transverse (x-y) 
plane at impact parameter b we use the Glauber Model 
density 

, n= (1 - o)pwn{x, y, b) + 2ap B (x, y, b) 
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where pwn and ps are the densities of wounded nu- 
cleon and binary collisions respectively, a = 0.145. 
The optical Glauber Model densities are obtained with 
Wood-Saxon densities for the Au nuclei pws( r ) — 
p /(exp((r-i? o )/a) + l) (p = 0.169f m - 3 , R a = 
6.38fm, a = 0.535fm) and the inelastic cross section is 
42mb. The energy density at the center of the fire- 
ball eo for b = is adjusted to reproduce the parti- 
cle multiplicity in the most central (0-5%) collisions in 
ideal hydrodynamic simulations. The initial density for 
other centralities is taken from the formula (3.5) with- 
out changing any parameters. In viscous hydrodynamic 
calculations the initial density is rescaled to take into 
account the additional entropy produced. In Fig. 2 is 
shown the entropy production in the different hydrody- 
namic evolutions. The entropy is normalized to the en- 
tropy in the ideal fluid simulation. In viscous hydro- 
dynamics the entropy increases with time, we chose to 
normalize the entropy in all the calculations to the same 
value at r — To = 6fm/c. This procedure yields, after 
hadronization, similar particle multiplicities in all the 
calculations. Entropy is produced mainly in the QGP 
phase AS/S ~ 20%, whereas in the hadronic matter its 
relative increase is only 2-3%. 

The freeze-out temperature is fixed to reproduce the 
transverse momentum spectra of pions in central colli- 
sions. The lifetime of the fireball is determined by the 
initial temperature, the expansion rate and the freeze- 
out temperature. The interplay of those effects makes 
the lifetime in the ideal fluid and vQGP+vHG scenarios 
very similar. The detailed shape of the freeze-out hyper- 
surface depends however on the amount of the accumu- 
lated transverse flow (Fig. 3). This has consequences on 
the resulting HBT radii and in particular on the ratio 

Rout I Rside • 



IV. RESULTS 

Transverse momentum spectra of pions are shown in 
Fig. 4. The freeze-out temperature is adjusted for each 
of the considered scenarios to reproduce pion spectra 
in the most central collisions for p± < 1.2GeV/c. In 
the ideal fluid expansion, reducing the freeze-out tem- 
perature means that the fluid expands longer and more 
transverse flow builds up. This effect dominates over 
the reduction of the final temperature and the spec- 
tra become harder. For the chosen initial conditions, 
Tp = 140MeV is optimal for the ideal fluid expanding 
from To = 0.6fm/c. Shear viscosity corrections in the 
plasma phase (scenario vQGP) result in additional trans- 
verse pressure in the early stage of the expansion. To 
reproduce the observed pion spectra the evolution must 
be shortened giving Tp — 150MeV. The situation is very 
different if dissipative corrections in the hadronic phase 
are allowed for (scenarios vHG or vQGH+vHG). Bulk 
viscosity leads to a substantial softening of light particle 
spectra, hydrodynamic evolution must be effective for a 



10 3 



> 
a> 
O 

>; 

73 

■a 

1- 
Q. 

CM 



10 -3U 



10"' 



10 ! 





7l + 


^\ ^\0-5% 










♦ 


40-50%\ S 

(x10- 6 )X. 

id. fl. 

- vQGP 

vHG 

vQGP+vHG 


- . . 



0.5 1 



1.5 



2 2.5 3 
p r [GeV/c] 



FIG. 4: (Color online) n + transverse momentum spectra spec- 
tra for Au-Au collisions at y^s = 200GeV and centralities 
0-5%, 5-10%, 10-15%, 15-20%, 20-30%, 30-40% and 40-50% 
(successively scaled down by powers of 0.1). The solid, dot- 
ted, dashed and dash-dotted lines represent the results of ideal 
hydrodynamic, vQGP, vHG and vQGP+vHG calculations re- 
spectively. Data are from PHENIX Collaboration [56]. 




FIG. 5: (Color online) Same as Fig 4 but for if -1 
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FIG. 6: (Color online) Same as Fig 4 but for protons. 



longer time in order to reproduce the p± spectra of pions. 
Depending on the amount of collective transverse flow ac- 
cumulated in the early phase of the dynamics it results 
in freeze-out temperatures 130-135MeV. Bulk viscosity 
corrections (Eq. 2.11) grow with the momentum of the 
particle and eventually become as large as the equilib- 
rium distribution / , it means that the formalism breaks 
down. Using the average bulk viscosity corrections at 
the freeze-out hypersurface we estimate that viscous hy- 
drodynamics with statistical emission of particles breaks 
down for pion momenta of 1.5GeV/c in the fluid rest 
frame. Pion spectra at large transverse momenta cannot 
be reliably described in the formalism used in this work. 
After adjusting the freeze-out conditions to reproduce 
pion spectra at soft momenta in central collisions, all ob- 
servables at different centralities are calculated without 
modifying the parameters of the model. We observe that 
pion spectra at different centralities are well described 
for p ± < 1.2GeV/c. 

In Figs. 5 and 6 we show the spectra of K + and pro- 
tons at different centralities. A first observation is that 
the slopes of the spectra for heavier particles obtained in 
scenarios with or without bulk viscosity do not differ as 
much as for pions. It is a consequence of the mass de- 
pendence of the bulk viscosity corrections in Eq. (2.11). 
Kaon production is overpredicted by hydrodynamic cal- 
culations in peripheral collisions which may be a manifes- 
tation of partial equilibration of strangeness [58] or of a 
nontrivial dependence of the thermal source size on cen- 
trality [59, 60]. The effective slopes of proton spectra for 
p± < 2GeV/c are well reproduced by all the calculations. 
The multiplicity of protons, reflected in the normaliza- 
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FIG. 7: (Color online) Elliptic flow coefficient for 7r + and K + 
(triangles and solid lines) and for protons (circles and dashed 
lines). Thin lines are for the ideal fluid model (all panels) 
and the thick lines for the various viscous hydrodynamic cal- 
culations, data are from the PHENIX Collab. [57]. In the 
upper, middle and lower panel are shown results obtained 
from vQGP+vHG, vHG and vQGP scenarios respectively. 



tion of the spectra in Fig. 6, is better described if bulk 
viscosity is present. The chemical freeze-out tempera- 
ture fitted from the particle number ratios is 165MeV 
[61, 62], significantly larger than the freeze-out temper- 
atures we use. Simulations where particles are emitted 
without bulk viscosity corrections (id. fl. or vQGP) un- 
derpredict the number of heavier particles. Bulk viscosity 
corrections reduce the number of light particles and lead 
to an increase in heavy particle production, resulting in 
an effective chemical non-equilibrium at freeze-out. Con- 
sequently, simulations including moderate bulk viscosity 
in the hadronic stage reproduce the proton number in 
spite of lower freeze-out temperatures. 

An important characteristic of the dynamics and of 
the equation of state of the fireball is the elliptic flow 
coefficient [6] . Most of the elliptic flow is created in the 
early phase of the expansion, and so the flow probes pres- 
sure gradients at that time. However at densities where 
freeze-out occurs the elliptic flow is still increasing dur- 
ing the hydrodynamic evolution. Extracting the shear 
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viscosity from the comparison of model calculations to 
the data requires a very precise, independent determina- 
tion of the freeze-out time. A strong constraint on the 
final density at freeze-out is given by the transverse mo- 
mentum dependence of vi for different species [32] and in 
particular the difference in the flow of pions and protons. 
This picture is more complicated if viscosity corrections 
at freeze-out are important. First, non-equilibrium cor- 
rections from shear viscosity could in principle be very 
different for different particles [55] and second, bulk vis- 
cosity corrections (Eq. 2.11) depend on the particle mass. 

In Fig. 7 we show the momentum dependent elliptic 
flow coefficient for light mesons and protons. The ideal 
fluid simulation (thin lines in all the panels) does not re- 
produce the meson-proton splitting present in the data. 
The elliptic flow of protons is too large. The same is 
true for the scenario where the viscosity is negligible in 
the hadronic phase (lower panel). Nonequilibrium cor- 
rections at freeze-out (both from shear and bulk viscosi- 
ties) are essentially zero in that case. Shear viscosity 
in the plasma phase changes the flow pattern reducing 
velocity gradients and leading to a decrease of the final 
elliptic flow. Meson and proton elliptic flow gets reduced 
in a similar way by the shear viscosity, and we cannot 
get enough meson-proton splitting. The situation is very 
different if we allow for viscosity corrections in the final 
stage of the expansion. The most important difference 
comes from corrections to the distribution functions at 
freeze-out. Shear viscosity corrections lead to an addi- 
tional reduction of the elliptic flow. However, bulk vis- 
cosity corrections reduce the transverse momenta of light 
mesons and the differential elliptic flow in p± is increased 
(two upper panels in Fig. 7). The same effect has been 
noticed in the estimates of bulk viscosity corrections at 
freeze-out in Ref. [53]. Bulk viscosity corrections are 
much smaller for heavy particles and are not sufficient to 
increase the value of the elliptic flow for protons. This 
and the lower freeze-out temperatures in the scenarios 
with hadronic bulk viscosity bring the species dependent 
elliptic flow to an agreement with the data. 

An interesting experimental observation is the mass 
scaling of identified particle elliptic flow in transverse 
mass [63] at small momenta. The mass ordering of the 
elliptic flow indicates a hydrodynamic origin of the ob- 
served flow. In Fig. 8 we present the elliptic flow as 
function of transverse mass for several identified hadrons. 
The results for different centralities are scaled by the ini- 
tial eccentricity of the fireball. The results in different 
panels correspond to different scenarios for the viscosi- 
ties. The best agreement with the transverse mass scaling 
is seen in the calculations with small freeze-out tempera- 
tures, vHG or vQGP+vHG. Our calculations in the ideal 
fluid or vQGP scenarios cannot reproduce the observed 
mass ordering. It must be noted that ideal fluid simula- 
tions with low freeze-out temperatures capture correctly 
the hydrodynamic origin of the mass ordering of the flow 
[65] . Similar results are obtained in hydrodynamics with 
shear viscosity for one value of r)/s [66] . 
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FIG. 8: (Color online) Elliptic flow as function of transverse 
mass from hydrodynamic calculation (scaled by the initial ec- 
centricity in the calculation, lines) and observed experimen- 
tally (scaled by the participant eccentricity, symbols) for 
(triangles and dotted line), A+A (circles and short dashed 
line), S+S (squares and solid line), all at centralities 10 — 40%, 
STAR Collab. data [63], for tt+ at centralities 10 - 20% (di- 
amonds and long dashed line), PHENIX Collab. data [57], 
and for protons at centralities 20 — 40% (reversed triangles 
and dash-dotted line), STAR Collab. data [64]. 
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FIG. 9: (Color online) Elliptic flow coefficient for charged par- 
ticles for different centralities. Stars represent ideal fluid re- 
sults, reversed triangles, triangles and circles represent the re- 
sults of vQGP, vHG and vQGP+vHG viscous hydrodynamic 
calculations respectively. The solid line with squares denotes 
the results of calculations using t/hg/s = 0.24, t/qgp/s = 0.08 
and the same bulk viscosity as in the scenarios vHG and 
vQGP+vHG. Data are from the STAR Collab. [64]. 
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FIG. 10: (Color online) Elliptic flow of mesons (triangles) and 
protons (circles) as function of transverse momentum [57], 
compared to model calculations. We present three simula- 
tions assuming an increasing shear viscosity to entropy den- 
sity ratio with decreasing temperature (vHG, tiqgp/s = 0, 
Vhg/s = 0.1, dashed lines), a constant one (t)qgp/s = 0.1, 
tihg/s — 0.1, dashed-dotted lines) and a decreasing one 
(vQGP+vHG, tiqgp/s = 0.16, r/Hc/s = 0.1, solid lines). 
Also are shown results for a calculation with a minimal QGP 
viscosity t/qgp/s = 0.08 and moderately large hadronic vis- 
cosity tjhg/s = 0.24 (dotted lines). Thick and thin lines 
represent meson and proton elliptic flow respectively. 



In Fig. 9 is plotted the average elliptic flow coefficient 
of charged particles at different centralities. The ideal 
fluid calculation overpredicts the elliptic flow in periph- 
eral collisions. The discrepancy increases with the im- 
pact parameter indicating that corrections to the ideal 
fluid dynamics should be more important in collisions 
where the hadronic phase is relatively more important. 
Comparing the three calculations with viscosities, we find 
that adding more dissipative mechanism reduces the fi- 
nal elliptic flow. However, the differences between the 
two scenarios with or without viscosity in the plasma 
phase are not very big. Most of the effect of the re- 
duction of the azimuthal asymmetry of the flow comes 
from the hadronic dissipation. One must conclude that 
the sensitivity of the elliptic flow to the shear viscos- 
ity in the early phase is strongly reduced if additional 
dissipation occurs below T c . It must be stressed that 
the assumed strength of shear viscosity in the hadronic 
medium is small. We performed also simulations in a sce- 
nario with a larger value of the viscosity in the hadron 
phase t]hg/ s = 0.24, with t]qgp/s = 0.08 and a freeze- 
out temperature of 135MeV. The result is plotted in Fig. 
9 (solid line with squares), the elliptic flow of charged 
particles is below the experimental values, which means 
that the assumed shear viscosity is too large. This is in 
line with previous calculations using a constant value of 
shear viscosity, where the simulations with rj/s = 0.08 
best reproduce the data [67]. 

In Fig. 10 is presented a direct comparison of the ellip- 
tic flow coefficient as function of transverse momentum 
for different choices of the temperature dependence of the 
shear viscosity to entropy density ratio. We take three 
different values of the QGP viscosity tjqgp/s — 0, 0.1 
and 0.16 and t]hgI s = 0.1. It means that we check 
three qualitatively different scenarios with rj/s increas- 
ing, constant or decreasing when the temperature drops 
below T c . The first observation is that the differences 
between the three calculations are small, moreover in- 
creasing rjQGp/s always leads to a decrease of the elliptic 
flow. We find a satisfactory description of the data with 
a small value of tjhg/s — 0.1. For a calculation using 
a minimal QGP shear viscosity tjqgp/s — 0.08 and a 
larger value of t/hg = 0.24 the calculated final elliptic 
flow of mesons is significantly below the data. Increasing 
VQGp/s leads to a decrease of V2, also increasing ijhg/s 
from 0.1 to 0.24 gives a strong reduction of the elliptic 
flow. Comparing these results to experimental data we 
obtain the following conclusions : both viscosity effects in 
the plasma and in the hadronic phase lead to a decrease of 
i>2, when choosing a small value oir/qcp/s = 0.08-0.1 the 
best results are obtained for a small value of t\hg! s = 0.1 
(and not 0.24), even if the QGP viscosity is zero the pre- 
ferred value of the hadronic viscosity is small. 

We calculate HBT correlation radii of pions emitted 
from the fireball in the most central collisions. The two 
particle correlation function sums all pairs of identical 
pions with interference effects [69, 70]. For a given total 
momentum of the pair k± the three-dimensional correla- 
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FIG. 11: (Color online) HBT radii for Au-Au collisions at 
centrality — 5%. Ideal fluid calculation (solid lines), viscous 
hydrodynamic models vQGP+vHG (dash-dotted lines), vHG 
(dashed lines), vQGP (dotted lines) and STAR Collab. data 
[68] (squares) are shown. 



tion function in the pion relative momentum is fitted with 
the Bertsch-Pratt formula [71, 72]. All the scenarios of 
the hydrodynamic expansion studied lead to HBT radii 
that are quite close to the data (Fig. 11). It is a con- 
sequence of the hard equation of state used, with only 
a minimal softening around T c [11]. The scenarios dif- 
fer by the freeze-out temperatures, a higher freeze-out 
temperature means a shorter lifetime and hence smaller 
values of the radii. R S ide measuring the geometrical size 
of the system at freeze-out decreases monotonically with 
increasing Tp. The description of the experimentally ob- 
served small value of the ratio Rout / Rside requires the 
use of a hard equation of state, early initial time of the 
expansion, dissipative effects and/or a Gaussian initial 
profile [11, 73]. Ideal fluid expansion and the expansion 
with only hadronic dissipation (solid and dashed lines in 
Fig. 11) give similar results for R ou t/ Rside- The same 
is true for scenarios with the same viscosity in the early 
stage but different hadronic dissipation (dotted and dash- 
dotted lines overlap in the lowest panel in Fig. 11). We 
can conclude that the ratio R ou t / Rside is sensitive to the 
early build up of the transverse flow and is not sensitive 
to viscosity effects at freeze-out. This is in contrast to 
the elliptic flow which is sensitive to dissipative effects 
at all the stages of the expansion. We cannot reproduce 
exactly the observed HBT radii, this may indicate that 



the amount of the early transverse flow in the expansion 
is too low. 



V. CONCLUSIONS 

We present a study of viscosity effects at different 
stages of the expansion of the fireball created in rela- 
tivistic heavy ion collisions. We introduce the possibility 
to have two different shear viscosities in the QGP and 
hadronic phases of the matter. Assuming zero shear vis- 
cosity or rj/s = 0.16 in the plasma, we test its impact on 
the final observables. The sensitivity of the final ellip- 
tic flow observables to the early viscosity is reduced by 
the dissipative effects in the hadronic phase. A crucial 
effect is the introduction of a moderate value of the bulk 
viscosity in the hadronic medium. Such an assumption 
is natural in a system with partial equilibration and fi- 
nite particle masses. The bulk and shear viscosities are 
treated as free parameters. We use however a relaxation 
time approximation to relate the values of the shear and 
bulk viscosities in the hadron gas. A moderate value of 
the bulk viscosity Q/s = 0.03-0.04 corresponds to a rela- 
tively small value of the shear viscosity rj/s = 0.1. Bulk 
viscosity in the late stage leads to a shift of the freeze-out 
temperature to a value allowing for a satisfactory descrip- 
tion of the mass ordering of the elliptic flow of identified 
hadrons without spoiling the agreement in the HBT radii 
and transverse momentum spectra. 

From the simulations using Glauber Model initial den- 
sities here presented we can conclude that the shear 
viscosity in the hadronic phase is in the range 0.1 < 
Vhg/ s < 0.24. Shear viscosity in the plasma phase leads 
to a decrease of elliptic flow, but it is the value of 77 at 
late stages that is the most important for the suppres- 
sion of the elliptic flow in all the cases vqgp/s < vhg/s, 
Vqgp I s = Vhg/s, or riQ G p/s > rj H G/s. Even taking 
Vqgp — leads to a preferred value of vhg/ s = 0.1 that 
best reproduces the data. Therefore, our results point- 
ing to a small value of the shear viscosity in the hadronic 
phase, are consistent with previous calculations using the 
same small value of r\ in the plasma and in the hadron 
fluids. 

The extracted shear viscosity is significantly below 
the microscopic estimates of the shear viscosity in the 
hadronic matter vhg/s — 1. Using rj/s ~ 1 in hydro- 
dynamic simulations with the large velocity gradients in 
heavy-ion collisions is beyond the range of applicability of 
Israel- Steward formalism, it would lead to severe numeri- 
cal instabilities, and is disfavored by elliptic flow data. It 
is not clear to the author what is the mechanism that pre- 
vents the expected large shear viscosity in the hadronic 
matter to become effective in the viscous hydrodynamic 
evolution of heavy ion collisions. Before concluding that 
the hadronic matter is indeed a low viscosity fluid, it 
should be checked whether the difference between trans- 
port model estimates and hydrodynamics is not due to a 
deficiency of the sudden freeze-out procedure used in the 
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calculations [74] or to a strong increase of the relaxation 
time t„ below T c analogously to what has been discussed 
for bulk viscosity in Ref. [50]. 

We find that the elliptic flow coefficient is significantly 
reduced due to viscosity effects both in the plasma and 
in the hadronic matter. It means that the extraction of 
the shear viscosity in QGP is difficult and can be reliably 
addressed only after precisely constraining the freeze-out 
conditions. By this we mean determining both the freeze- 
out temperature and the nonequilibrium effects in the fi- 
nal state. It is interesting to note that the HBT radii have 
a simple dependence on the choice of the freeze-out. The 
radii increase for a larger lifetime of the system, caused by 
a smaller freeze-out temperature. The ratio Rout/Rside 
is almost insensitive to the freeze-out condition, but it 
depends on the amount of the transverse flow generated 
in the early phase. One of the mechanism increasing the 



early transverse flow is due to the shear viscosity in the 
plasma phase. 

Let us close by repeating the observation that the 
introduction of bulk viscosity in the hadronic medium 
changes the freeze-out conditions in the hydrodynamic 
expansion of the fireball. This allows for a good and si- 
multaneous description of transverse momentum spectra, 
identified particle elliptic flow and HBT radii. 
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